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4  Introduction 

The  early  detection  of  breast  cancer  reduces  patient  mortality.  The  most  common  method  of 
breast  cancer  detection  is  palpation.  However,  lesions  that  lie  deep  within  the  breast  are  dif¬ 
ficult  to  palpate  when  they  are  small.  Thus,  a  method  of  Remote  Palpation,  which  may  allow 
the  detection  of  small  lesions  lying  deep  within  the  breast,  is  currently  under  investigation. 
In  this  method,  acoustic  radiation  force  is  used  to  apply  localized  (on  the  order  of  2  mm3) 
forces  within  tissue,  and  the  resulting  tissue  displacements  are  mapped  using  ultrasonic  cor¬ 
relation  based  methods.  An  area  that  is  stiffer  than  the  surrounding  medium  distributes  the 
force,  resulting  in  larger  regions  of  displacement,  and  smaller  maximum  displacements.  The 
resulting  displacement  maps  are  directly  correlated  with  local  tissue  stiffness.  In  the  second 
year  of  support,  we  have  expanded  our  finite  element  models  of  tissue  mechanical  response 
to  radiation  force;  we  have  developed  finite  element  models  of  tissue  thermal  response  to 
acoustic  radiation  force;  we  have  performed  experiments  that  validate  the  mechanical  finite 
element  model;  we  have  constructed  calibrated  tissue  phantoms  and  methods  for  measuring 
their  mechanical  properties,  we  have  implemented  acoustic  Remote  Palpation  using  multiple 
pushing  locations  in  phantoms,  and  we  have  obtained  initial  clinical  data. 


5  Body 

We  have  developed  a  finite  element  model  of  Remote  Palpation,  which  has  been  validated 
experimentally  [2,  1,  3].  The  results  of  our  phantom  experiments  and  simulations  are  sum¬ 
marized  in  Appendices  10.1  and  10.2.  This  represents  completion  of  Task  1  of  the  statement 
of  work. 

We  have  fabricated  multiple  phantoms  with  varying  Young’s  modulii,  with  embedded  lesions 
of  varying  sizes,  and  we  have  measured  the  Young’s  moduli  of  these  phantoms.  We  have 
developed  pulse  sequences  on  the  Siemens  Elegra  scanner  to  implement  Remote  Palpation  in 
a  single  pushing  location,  and  we  have  utilized  an  automated  translation  stage  to  translate 
the  transducer  in  order  to  use  this  pulse  sequence  to  interrogate  an  extended  Field-of-View 
(FOV).  Using  the  experimental  setup,  we  have  implemented  Remote  Palpation  in  phantoms 
using  multiple  pushing  locations,  and  created  the  first  Remote  Palpation  image  of  an  em¬ 
bedded  lesion.  This  image  exhibited  increased  contrast  over  the  corresponding  ultrasonic 
B-mode  image  [3]  (Appendix  10.2).  These  experiments  and  results  represent  completion  of 
Task  2  of  the  statement  of  work. 

We  have  developed  a  finite  element  model  of  the  thermal  response  of  tissue  to  acoustic  ra¬ 
diation  force,  and  have  performed  initial  simulations  (Figure  1).  These  simulations  indicate 
that  the  spatial  extent  of  the  heating  is  limited  to  the  focal  region  of  the  acoustic  beam,  and 
that  the  maximum  increases  in  temperature  will  be  less  than  1°C  during  method  implemen¬ 
tation  1.  This  represents  partial  completion  of  Task  3  of  the  statement  of  work.  Further 
simulations  are  required  to  investigate  the  different  pulsing  sequences  that  will  be  developed 
for  further  clinical  implementation. 

We  have  initiated  clinical  trials  of  the  Remote  Palpation  method.  Three  patients  have  been 
imaged  to  date.  Detectable  tissue  displacements  were  generated  in  one  of  the  patients  (~  6 
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Figure  1:  Top:  Simulated  displacement  profile,  and  Bottom:  simulated  temperature  increases 
resulting  from  a  typical  remote  palpation  beam  sequence  (in  vitro  Iapta  =  250  W/cm2,  Tissue 
Young’s  modulus  =  3.0  kPa).  Note  the  different  spatial  scales  on  the  two  images.  The  yellow 
contour  in  the  temperature  plot  is  approximately  the  size  of  the  pushing  region  of  the  acoustic 
beam.  The  shape  of  the  white  region  in  the  displacement  plot  is  similar  to  that  of  the  focal  region, 
although  it  is  more  broad.  This  is  due  to  the  shear  stresses  present  in  the  medium  ( i.e .  the  tissue 
’drags’  the  regions  of  tissue  that  surround  the  region  where  the  force  is  applied).  The  relatively 
small  temperature  increases  (light  blue  through  red  regions,  maximum  increase  0.43°  C)  are  limited 
to  the  focal  region  of  the  acoustic  beam.  The  transducer  is  on  the  left  side  of  these  images. 

microns,  Figure  2).  This  comprises  partial  completion  of  Task  4  of  the  statement  of  work. 
We  have  determined  from  these  trials  that  further  modifications  are  required  in  the  beam 
sequences  and  the  system  power  supplies  in  order  to  generate  the  required  acoustic  radiation 
forces  in  breast  tissue.  Even  so,  we  are  encouraged  by  the  results  to  date,  and  are  working 
to  overcome  the  challenges  that  have  been  identified. 

6  Key  Research  Accomplishments 

•  Enhanced  finite  element  model  of  radiation  force  induced  tissue  motion. 

•  Developed  finite  element  model  of  the  thermal  response  of  tissue  to  radiation  force. 

•  Developed  calibrated  tissue  phantoms. 

•  Developed  measurement  methods  and  measured  the  mechanical  properties  of  the  tissue 
phantoms. 

•  Developed  experimental  pulse  sequences  to  implement  Remote  Palpation  using  multi¬ 
ple  pushing  locations. 

•  Implemented  Remote  Palpation  with  multiple  pushing  locations  in  phantom. 
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lateral  position  (mm) 
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Figure  2:  TOP:  In  vivo  displacement  image  after  10  milliseconds  of  pushing  in  a  single  pushing 
location  in:  a)  breast  (33  year  old  female).  The  tracking  beams  were  spaced  0.6  mm  apart,  and 
were  sequenced  to  track  laterally  in  the  center,  then  out  to  the  right,  then  back  from  center  to 
left.  The  tissue  begins  to  relax  during  the  tracking  sequences,  which  is  the  reason  for  the  larger 
displacements  on  the  right  side  than  the  left.  BOTTOM:  Corresponding  displacement  through  time 
at  the  focal  point  of  the  transducer.  This  image  constitutes  the  first  in  vivo  RP  image.  Note 
that  this  displacement  was  achieved  at  a  pushing  pulse  duty  cycle  of  3.4%.  Significantly  higher 
displacements  are  expected  for  an  equivalent  power,  100%  duty  cycle  pushing  pulse.  The  estimated 
intensity  of  the  pushing  beams  was  240  W/cm2,  and  an  F/l  focal  configuration  was  used.  These 
experiments  were  performed  using  a  non-attenuating  standoff  pad  (10  mm). 


•  Generated  the  first  Remote  Palpation  images  using  this  data. 

•  Initiated  clinical  trials  of  the  method. 


7  Reportable  Outcomes 

•  Paper  published  in  Ultrasonic  Imaging,  entitled:  “A  finite  element  model  of  Remote 
Palpation  of  breast  lesions  using  radiation  force:  factors  affecting  tissue  displacement.” 

•  Paper  submission  to  Journal  of  the  Acoustical  Society  of  America,  entitled:  “On  the 
feasibility  of  remote  palpation  using  acoustic  radiation  force” . 

•  Poster  at  the  DOD  BCRP  Era  of  Hope  meeting,  in  Atlanta,  GA  (June  8-11,  2000) 
entitled:  “Acoustic  Remote  Palpation:  Initial  Simulation  and  Experimental  Results” 

•  Presentation  at  the  Twenty-fifth  international  symposium  on  ultrasonic  imaging  and 
tissue  characterization  in  Arlington,  Va  (May  24,  2000)  entitled:  “An  experimental 
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investigation  of  the  Required  Acoustic  Power  for  In  Vivo  Implementation  of  Radiation 
Force  Based  Imaging.” 


•  Presentation  at  the  2000  Annual  Convention  of  the  American  Institute  of  Ultrasound 
in  Medicine  in  San  Francisco,  CA  (April  5,  2000)  entitled:  “Acoustic  radiation  force 
induced  motion  to  characterize  variations  in  tissue  stiffness:  initial  experimental  re¬ 
sults.” 

•  Presentation  at  the  1999  IEEE  International  Ultrasonics  symposium  at  Lake  Tahoe, 
NV  (October  19,  1999)  entitled:  “Finite  element  analysis  of  radiation  force  induced 
tissue  motion  with  experimental  validation.” 

•  Patent  application  filed,  entitled:  “Method  and  Apparatus  for  the  identification  and 
characterization  of  regions  of  altered  stiffness” . 

•  Whitaker  grant  application  filed,  and  awarded,  with  funding  to  begin  in  September 
2001,  entitled:  “A  Radiation  Force  Based  Ultrasonic  Imaging  System  for  the  Early 
Detection  of  Breast  Cancer” . 


8  Conclusions 

The  clinical  implications  of  the  results  obtained  during  the  second  year  of  support  are  promis¬ 
ing.  The  simulations  and  experimental  results  are  consistent,  and  provide  the  framework  for 
the  design  of  a  Remote  Palpation  imaging  system.  The  first  Remote  Palpation  image  was 
created  using  multiple  pushing  locations  to  interrogate  an  extended  field-of-view  (FOV)  in  a 
calibrated  phantom  with  and  embedded  lesion.  This  image  showed  increased  contrast  over 
the  matched  ultrasonic  B-mode  image.  The  initial  thermal  model  results  indicate  that  in¬ 
creases  in  tissue  temperature  during  method  implementation  will  be  less  than  1°C,  which  is 
within  the  FDA  guidelines.  The  initial  clinical  data  indicates  that  increases  in  acoustic  power 
output  and  modifications  in  the  temporal  power  distribution  will  be  required  for  clinical  im¬ 
plementation  of  the  method.  Key  issues  requiring  further  investigation  include:  thermal  and 
power  issues,  development  of  pulse  sequences  with  varying  temporal  power  distributions,  and 
investigation  of  methods  of  information  display. 


7 


9  Bibliography 
References 


[1]  K.R.  Nightingale,  R.W.  Nightingale,  M.L.  Palmeri,  and  G.E.  Trahey.  Finite  element  anal¬ 
ysis  of  radiation  force  induced  tissue  motion  with  experimental  validation.  In  Proceedings 
of  the  1999  IEEE  Ultrasonics  Symposium ,  pages  1319-1323,  1999. 

[2]  K.R.  Nightingale,  R.W.  Nightingale,  M.L.  Palmeri,  and  G.E.  Trahey.  A  finite  element 
model  of  remote  palpation  of  breast  lesions  using  radiation  force:  Factors  affecting  tissue 
displacement.  Ultrasonic  Imaging ,  22(1)  :35— 54,  January  2000. 

[3]  K.R.  Nightingale,  M.L.  Palmeri,  R.W.  Nightingale,  and  G.E.  Trahey.  On  the  feasibility 
of  remote  palpation  using  acoustic  radiation  force.  J.  Acoust.  Soc.  Am.,  in  submission. 


8 


10  Appendices 


10.1  Proceedings  from  IEEE  UFFC  meeting 


10 


FINITE  ELEMENT  ANALYSIS  OF  RADIATION  FORCE  INDUCED  TISSUE  MOTION 

WITH  EXPERIMENTAL  VALIDATION 


Kathryn  Nightingale,  Roger  Nightingale,  Mark  Palmeri,  and  Gregg  Trahey 


Duke  University  Department  of  Biomedical  Engineering 
Durham,  NC  27708 
kathy.nightingale@duke.edu 


ABSTRACT 

An  ultrasonic  radiation  force-based  method  for  remote  palpation 
of  tissue  is  investigated.  The  use  of  radiation  force  to  image  tis¬ 
sue  stiffness  has  been  proposed  by  several  researchers.  In  this 
paper,  the  potential  for  using  a  diagnostic  ultrasound  system  to 
both  apply  radiation  force  and  track  the  resulting  tissue  displace¬ 
ments  is  investigated  using  Finite  Element  Methods  (FEM),  and 
the  results  are  compared  with  experimental  results.  Remote  pal¬ 
pation  is  accomplished  by  interspersing  high  intensity  pushing 
beams  with  low  intensity  tracking  beams.  This  generates  local¬ 
ized  radiation  forces  which  can  be  applied  throughout  the  tissue, 
with  the  resulting  displacement  patterns  determined  using  cor¬ 
relation  techniques.  An  area  that  is  stiffer  than  the  surrounding 
medium  distributes  the  force,  resulting  in  larger  regions  of  dis¬ 
placement.  and  smaller  maximum  displacements.  The  resulting 
displacement  maps  provide  information  as  to  the  location  and  size 
of  regions  of  increased  stiffness. 

We  have  developed  an  FEM  model  that  predicts  displace¬ 
ments  resulting  from  acoustic  radiation  force  fields  generated  by 
diagnostic  transducers  in  various  complex  media.  We  perform 
a  parametric  analysis  of  varying  tissue  and  acoustic  beam  char¬ 
acteristics  on  radiation  force  induced  tissue  displacements.  Dis¬ 
placements  are  on  the  order  of  microns,  with  considerable  dif¬ 
ferences  in  displacement  patterns  in  the  presence  and  absence  of 
a  lesion  (or  stiff  inclusion).  Initial  experimental  results  are  pre¬ 
sented  that  support  the  findings  in  the  model. 

1.  INTRODUCTION 

One  of  the  most  successful  methods  of  early  breast  cancer  detec¬ 
tion  is  palpation.  This  is  because  the  Young's  modulus  (Le.  stiff¬ 
ness)  of  malignant  tumors  is  often  an  order  of  magnitude  greater 
than  that  of  normal  breast  tissue  [15,  5].  Therefore,  cancerous  le¬ 
sions  feel  ‘harder’  or  ‘stony’  as  compared  to  normal  breast  tissue. 
Several  imaging  modalities  are  under  investigation  that  provide 
information  about  the  stiffness  of  tissue.  These  include  eiastog- 
raphv.  in  which  local  variations  in  tissue  strain  are  determined 
by  measuring  local  displacements  that  occur  during  global  tissue 
compression  [15,  12,  11],  and  sonoeiasticity,  in  which  low  fre¬ 
quency  shear  wave  propagation  is  imaged  using  Doppler  meth- 

This  work  was  supported  by  DOD/BCRP  grant  BC972755.  We  thank 
Dr.  Tim  Hall  of  the  University  of  Kansas  for  his  support  in  phantom 
design:  Dr.  Bruce  McDermott  and  Ms.  Cindy  Kerby  of  Siemens  Medical 
Systems.  Ultrasound  Group  for  their  system  support;  Altair  Computing, 
Inc.  for  their  meshing  and  post-processing  support,  and  Intel  Corporation 
for  their  technical  and  in-kind  support. 
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ods,  from  which  the  elastic  modulus  (Young’s  modulus)  of  tissue 
can  be  estimated  [22,  6,  13,  4,  7], 

The  use  of  acoustic  radiation  force  to  image  tissue  stiffness 
is  a  fairly  new  concept  that  appears  to  hold  some  promise  [14, 
20.  18,  9.  2].  All  of  the  proposed  methods  involve  using  radiation 
force  to  displace  tissue  in  a  remote  location,  and  monitoring  the 
tissue  response.  However,  they  differ  in  both  system  design  and 
tissue  response  characterization. 

Acoustic  radiation  force  is  a  phenomenon  associated  with  the 
propagation  of  acoustic  waves  through  a  dissipative  medium.  It  is 
caused  by  the  energy  density  gradient  that  occurs  in  the  medium, 
arising  either  from  absorption  or  reflection  of  the  wave  [19].  This 
gradient  results  in  the  application  of  a  force  in  the  direction  of 
wave  propagation.  In  an  absorbing  medium,  and  under  plane 
wave  assumptions,  this  force  can  be  represented  by  the  follow¬ 
ing  equation  [19.  10.  1.  17]: 

p .  If  absorbed  2.Q.I  ^  j  ^ 

C  C 

where  F  is  acoustic  radiation  force  [kg/(s2cm2)].  or  [dynes/(  1000 
cm3)].  Wabsorbed  [Watts/UOO  cm3)]  is  the  power  absorbed  by 
the  medium  at  a  given  point  in  space,  c  [m/s]  is  the  speed  of 
sound  in  the  medium,  a  [m-1]  is  the  absorption  coefficient  of  the 
medium,  and  I  [Watts/cm2]  is  the  temporal  average  intensity  at 
a  given  point  in  space.  Equation  1  provides  a  simple  relationship 
between  the  temporal  average  intensity  of  an  acoustic  beam,  and 
the  resulting  radiation  force  (which  is  in  the  form  of  a  body  force, 
or  force  per  unit  volume). 

The  method  under  investigation  herein  is  called  Remote  Pal¬ 
pation.  In  this  method,  a  single  transducer  on  a  diagnostic  ul¬ 
trasound  system  is  used  to  apply  localized  (on  the  order  of  2 
mm3)  radiation  forces  within  tissue,  and  to  track  the  resulting 
tissue  displacements  which  are  determined  using  ultrasonic  cor¬ 
relation  based  methods.  An  area  that  is  stiffer  than  the  surround¬ 
ing  medium  distributes  the  force,  resulting  in  larger  regions  of 
displacement,  and  smaller  maximum  displacements.  This  study 
was  designed  to  determine  the  feasibility  of  Remote  Palpation. 
A  Finite  Element  Model  (FEM)  is  presented  with  which  a  para¬ 
metric  analysis  of  varying  tissue  and  acoustic  beam  parameters 
is  performed.  Initial  experiments  designed  per  the  FEM  findings 
are  also  presented,  which  provide  good  agreement  with  the  sim¬ 
ulations. 

2.  SIMULATION  METHODS 

Model  implementation  is  performed  using  a  two  step  approach: 
first  the  spatially  distributed  intensity  field  from  a  given  trans- 
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Figure  1:  Left:  Finite  element  mesh.  The  transducer  is  on  the  left 
(back)  side  of  the  page,  with  its  axial  dimension  along  the  central 
axis  of  the  mesh.  The  diameter  of  the  outer  hemisphere  is  5  cm, 
and  is  intended  to  represent  an  effectively  infinite  boundary  con¬ 
dition.  Right:  Magnified  version  of  the  central  plane  of  the  mesh, 
showing  the  increased  element  density  in  the  focal  region,  the  ra¬ 
diation  force  field  contours  for  the  F/l  transducer  configuration, 
and  the  outline  of  a  spherical  1.0  cm  diameter  lesion. 


ducer  and  set  of  transmit  parameters  is  determined  and  the  asso¬ 
ciated  radiation  force  field  is  computed  using  Equation  1.  Second, 
finite  element  methods  are  used  to  solve  for  the  resulting  tissue 
displacement  patterns. 

FIELDJI  (http://www.it. dtu.dk/ — jaj/field/)  [3],  an  acoustic 
field  simulation  program,  was  used  to  model  the  temporal  aver¬ 
age  intensity  patterns  corresponding  to  various  diagnostic  trans¬ 
ducer  configurations.  The  transducer  specifications  were  selected 
to  model  the  Siemens  75L40  linear  transducer  array  (center  fre¬ 
quency  7.2  MHz,  element  height  5  mm.  aperture  width  dependent 
upon  the  number  of  active  elements:  Siemens  Vledical  Systems. 
Ultrasound  Group.  Issaquah,  WA).  The  simulated  voltage  used  to 
excite  each  element  was  fixed:  thus,  when  more  elements  were 
excited,  more  energy  was  transmitted.  The  transducer  was  fo¬ 
cused  laterally  at  the  location  of  its  elevation  focal  point  (which 
is  determined  by  its  lens). 

The  intensity  fields  for  each  transducer  configuration  (or  f- 
number.  which  is  the  parameter  we  use  as  an  indication  of  the 
number  of  elements  included  in  the  transmit  aperture)  were  com¬ 
puted  in  the  axial/lateral  plane  of  the  transducer  and  discretized 
into  three  contours  of  constant  intensity.  Once  the  contours  and 
the  desired  boundaries  were  in  place,  they  were  imported  into  a  fi¬ 
nite  element  mesh  generation  program  (Hvpermesh.  Altair  Com¬ 
puting  Inc.,  Troy,  MI).  The  mesh  was  generated  within  this  plane, 
and  then  rotated  90  degrees  around  its  central  axis,  in  order  to 
achieve  a  three  dimensional  model  (Fig.  1).  The  model  was  con¬ 
strained  on  the  surface  of  the  distal  quarter  of  the  hemisphere 
(opposite  the  transducer  location),  thus  modeling  a  breast  resting 
on  a  concave  platform.  Because  the  boundary  conditions  did  not 
allow  motion  normal  to  the  planar  surfaces  in  Figure  1,  symme¬ 
try  assumptions  hold,  and  solution  of  this  mesh  simulates  motion 
in  one  quarter  of  the  axi-symmetric  three  dimensional  spherical 
model.  For  the  parametric  analysis,  the  radiation  force  value  ap¬ 
plied  within  each  intensity  contour  was  computed  using  Eq.  (1). 
Peak  in  situ  intensities  of  90  W/cm2  and  15.3  W/cm2  were  used 
for  the  F/l  and  F/3  transducer  configurations,  respectively.  The 
speed  of  sound  was  1530  m/s,  and  the  absorption  coefficient  was 
.415  cm~\  or  0.5  dB/cm/MHz.  The  resulting  maximum  body 
forces  were  4882  and  830  dynes/cm3,  for  the  F/I  and  F/3  trans¬ 
ducer  configurations,  respectively.  The  density  was  1000  kg/m3. 
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Figure  2:  Contours  of  displacement  in  normal  tissue  (top. 
E=10kPa),  in  normal  tissue  with  a  0.5  cm  diameter  lesion  (mid¬ 
dle.  E=400kPa),  and  in  normal  tissue  with  a  1.0  cm  diameter  le¬ 
sion  (bottom.  E=400kPa).  The  transducer  is  located  on  the  left 
side  of  these  images,  and  the  profiles  are  symmetric  about  the 
x-axis.  The  location  of  the  lesions  is  highlighted  in  gray.  The  le¬ 
sions  are  moving  toward  the  right  as  rigid  bodies,  with  a  smaller 
maximum  displacement  associated  with  the  larger  lesion. 


and  the  Young’s  modulus  (E)  was  10  kPa.  When  lesions  of  vary¬ 
ing  stiffness  were  introduced,  the  Young's  moduli  ranged  from 
1  to  400  kPa.  which  is  within  the  potential  range  of  normal  and 
diseased  breast  tissue  [15,  5.  16].  The  tissue  was  assumed  to  be 
elastic. 

The  equations  of  motion  were  then  solved  numerically  using 
LS-DYNA  (Livermore  Software  Technology  Corporation.  Liver¬ 
more,  CA).  Radiation  force  was  applied  using  a  ramp  and  hold 
function,  where  the  ramp  up  to  the  maximum  force  value  was  lin¬ 
ear.  occurring  over  2  milliseconds,  and  the  force  was  then  applied 
continuously  for  40  milliseconds.  The  same  mesh  was  utilized  for 
all  simulations.  Forces  were  applied  to  different  mesh  elements 
depending  upon  the  transducer  configuration,  and  tissue  proper¬ 
ties  were  applied  to  different  elements  depending  upon  the  lesion 
size  being  modeled.  There  were  7440  elements  in  the  mesh.  The 
run  times  varied  with  Young’s  modulus,  and  were  anywhere  from 
5  minutes  to  6  hours  on  a  450  MHz  Pentium  PC. 

In  order  to  perform  experimental  validation  of  the  model,  it 
was  also  run  using  material  properties  corresponding  to  the  cali¬ 
bration  phantom  (a=0.53  dB/cm/MHz,  c=1550  m/s,  E=0.3  kPa), 
and  acoustic  parameters  corresponding  to  those  in  the  experi¬ 
ment. 
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Figure  3:  Magnified  view  of  contours  of  displacement  (dashed 
lines)  in  normal  breast  tissue  (E=IOkPa),  for  the  F/l  (top)  and  F/3 
transducer  configurations  (bottom).  The  transducer  is  on  the  left 
side  of  the  page.  The  outer  boundary  of  the  associated  radiation 
force  field  is  superimposed  on  each  figure  as  a  solid  black  line. 
Each  contour  level  represents  a  \09c  decrease  in  displacement. 
Note  that  the  shape  of  the  displacement  contours  are  similar  to 
the  shape  of  the  radiation  force  fields:  the  F/l  contours  are  shifted 
slightly  in  front  of  the  focus,  but  are  fairly  small,  whereas  the  F/3 
contours  are  much  longer  axially. 


3.  EXPERIMENTAL  METHODS 

Experiments  were  performed  with  a  Siemens  Elegra  scanner  (Sie¬ 
mens  Medical  Systems.  Ultrasound  Group.  Issaquah.  WA).  that 
has  been  modified  to  allow  user  control  of  the  acoustic  beam  se¬ 
quences  and  intensities,  as  well  as  providing  access  to  the  raw 
Radio-Frequency  (RF)  data.  The  beam  sequence  was  designed  to 
transmit  a  series  of  13  spatially  distributed  tracking  beams,  fol¬ 
lowed  by  10  milliseconds  of  high  intensity  ’pushing”  beams  fired 
along  a  single  line  of  flight  in  the  center  of  the  tracking  beams, 
followed  by  another  series  of  tracking  beams  interspersed  with 
pushing  beams  (every  other  beam).  This  sequence  was  repeated 
for  100  milliseconds,  and  the  raw  RF  data  was  stored  for  off-line 
processing.  Experiments  were  performed  in  both  an  F/l  (more 
elements)  and  F/3  (fewer  elements)  focal  configuration,  using  the 
same  transmit  parameters  in  both  cases.  Thus,  in  the  calibrated 
phantom  (F/3  configuration,  q-0.53  dB/cm/MHz.  c=  1550  m/s, 
and  Young's  modulus  =  0.3  kPa),  the  derated  spatial  peak  tem¬ 
poral  average  intensity  was  1.0  W/crrr,  and  the  Mechanical  In¬ 
dex  was  0.3.  In  the  lesion  phantom  (F/l  configuration,  a=0.32 
dB/cm/MHz.  c=l550  m/s.  Young's  modulus  =  0.05  kPa  (tissue) 
and  13.5  kPa  (lesion)),  the  derated  spatial  peak  temporal  average 
intensity  was  2.9  W/cm2,  and  the  Mechanical  Index  was  0.6. 

Off-line  data  processing  was  accomplished  by  performing 
cross-correlation  between  sequentially  acquired  tracking  lines. 
Each  tracking  line  was  divided  into  a  series  of  search  regions, 
and  the  location  of  the  peak  in  the  cross  correlation  function  be¬ 
tween  a  kernel  in  the  first  tracking  line  and  the  corresponding 
search  region  in  the  next  tracking  line  was  used  to  estimate  the 
tissue  displacement  in  that  region. 


Figure  4:  Simulated  results  for  variations  in  force  and  stiffness. 
These  portray  a  linear  relationship  between  force  and  peak  dis¬ 
placement.  and  an  inverse  relationship  between  stiffness  and  dis¬ 
placement. 


4.  RESULTS  AND  DISCUSSION 

Figure  2  provides  contours  of  constant  displacement  in  the  pres¬ 
ence  and  absence  of  lesions  of  varying  size  generated  using  an  F/l 
transducer  configuration.  Considerable  differences  in  displace¬ 
ment  profiles  exist  in  the  presence  and  absence  of  the  lesions.  The 
peak  displacement  in  the  absence  of  the  lesion  is  1.2  ^m.  whereas 
for  the  0.5  cm  lesion  the  peak  is  0.7/im,  and  for  the  1 .0  cm  lesion 
it  is  0.3/im.  In  addition,  the  lesions  appear  to  exhibit  rigid  body 
motion,  and  thus  the  lateral  extent  of  the  displacement  profiles  is 
much  larger  in  the  presence  of  a  lesion.  The  displacement  profile 
in  the  absence  of  a  lesion  is  much  narrower,  and  is  similar  to  the 
shape  of  the  radiation  force  field  (Fig.  3.  top).  Note  that  in  the  F/3 
case  (Fig.  3,  bottom),  the  displacement  profile  is  elongated  axi¬ 
ally,  which  results  from  the  larger  volume  of  the  radiation  force 
field. 

The  simulations  were  performed  for  varying  forces  and  tissue 
stiffnesses  (Figure  4).  As  expected,  a  linear  relationship  exists 
between  peak  displacement  and  force,  and  an  inverse  relationship 
is  observed  between  tissue  stiffness  and  peak  displacement. 

Excellent  agreement  is  evident  between  the  experimental  and 
simulation  results  (Fig.  5).  Both  simulation  and  experiment  por¬ 
tray  a  fairly  elongated  axial  displacement  pattern  in  the  F/3  trans¬ 
ducer  configuration,  with  comparable  lateral  extent,  resulting  in 
teardrop-shaped  profiles.  In  addition,  the  displacement  magni¬ 
tudes  are  similar  (maximum  displacements  of  1.4  jum  in  the  sim¬ 
ulation,  and  1.9  ^m  in  the  experiment,  after  20  milliseconds  of 
high  intensity  insonification). 

Experiments  were  also  performed  in  a  phantom  that  was  de¬ 
signed  with  an  embedded  0.8  cm  diameter  lesion.  This  lesion 
was  considerably  more  stiff  than  the  surrounding  tissue  (~  270 
times).  As  predicted  by  the  simulations,  the  displacement  maps 
are  different  in  the  presence  and  absence  of  the  lesion,  both  in 
magnitude  and  in  shape,  and  the  lesion  appears  to  be  moving  as 
a  rigid  body.  (Fig.  6).  Note  also  that  the  axial  extent  of  the  dis¬ 
placement  pattern  in  the  absence  of  a  lesion  is  much  smaller  in 
the  F/l  case  than  in  the  F/3  case,  as  predicted  by  the  simulations 
(Fig.  6,  left  and  Fig.  5,  right  vs.  Fig.  3) 
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Figure  5:  Simulated  (left)  and  experimental(right)  displacement 
maps  for  an  F/3  transducer  configuration  in  a  homogeneous  phan¬ 
tom  (E=0.3kPa).  after  20  milliseconds  of  high  intensity  insonifi- 
cation  {in  situ  Iap<a  =  l-0  W/crir,  Ml=0.3).  The  transducer  is  lo¬ 
cated  at  the  top  of  these  images.  The  shape  of  the  displacement 
profiles  is  similar,  although  the  location  of  the  peak  displacement 
is  slightly  closer  to  the  transducer  in  the  experimental  case. 


The  temporal  response  of  the  phantom  was  slightly  slower 
than  that  predicted  by  the  simulations  (Fig.  7),  which  may  be  due 
to  the  assumption  of  an  elastic  medium  in  the  simulations.  Of 
particular  interest  is  the  fact  that  in  the  presence  of  a  lesion,  the 
response  is  much  slower.  Thus,  after  even  5  milliseconds  of  force 
application  considerable  differences  in  displacement  profiles  ex¬ 
ist.  The  challenge  is  generating  enough  force  to  create  detectable 
displacements  {i.e.  greater  than  0.5  /i m  [21])  within  5  millisec¬ 
onds. 

The  question  of  the  maximum  intensity  output  from  diagnos¬ 
tic  scanners  is  complex,  Because  the  FDA  limits  are  considerably 
lower  than  the  maximum  capacity  of  current  diagnostic  scanners, 
manufacturers  do  not  quantify  this  parameter.  The  short  tempo¬ 
ral  requirements  are  also  an  issue  (only  5  milliseconds),  because 
most  system  measurements  are  performed  assuming  a  steady- 
state  application.  We  conservatively  estimate  that  for  a  clini¬ 
cally  relevant  range  of  tissue  stiffnesses,  the  required  temporal 
average  intensity  exceeds  that  currently  available  from  diagnos¬ 
tic  scanners  by  a  factor  of  30.  However,  when  one  considers  the 
increase  in  radiation  force  associated  with  nonlinear  propagation 
[14.  17,  8]  (we  have  observed  increases  in  radiation  force  by  a 
factor  of  2.6  in  breast  applications  in  vivo  [8]),  and  the  fact  that 
two-dimensional  transducers  with  more  surface  area  will  achieve 
higher  intensities  (simulations  suggest  a  factor  of  3),  this  factor 
drops  to  only  4.5. 

5.  CONCLUSIONS 

The  simulations  and  experiments  presented  herein  provide  the 
framework  for  the  design  of  a  Remote  Palpation  imaging  system. 
For  a  clinically  relevant  range  of  elastic  moduli  and  lesion  sizes, 
considerably  different  displacement  patterns  were  generated.  The 
forces  required  to  generate  these  displacements  are  higher  than 
those  currently  used  in  diagnostic  ultrasound.  They  should  not 
pose  a  thermal  danger  to  the  patient  due  to  their  short  duration, 
however,  they  may  require  the  design  of  more  efficient  transduc- 
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Figure  6:  Experimental  displacement  maps  after  20  millisec¬ 
onds  of  focusing  high  intensity  pulses  {in  situ  IJfp,<J=2.9  W/cnr, 
MI=0.6)  in  a  homogeneous  region  of  the  lesion  phantom  (left, 
E=0.05  kPa)  and  in  the  center  of  a  0.8  cm  diameter  lesion  in  the 
same  phantom  (right,  lesion  E=  1 3  kPa).  An  F/l  focal  config¬ 
uration  was  used,  with  the  transducer  located  at  the  top  of  the 
images. 
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Figure  7:  Simulated  (top)  and  experimental  (bottom)  normalized 
displacement  through  time  at  the  focal  point  of  the  transducer. 
The  results  are  normalized  to  the  displacement  after  40  millisec¬ 
onds  in  the  no  lesion  case.  The  temporal  response  of  the  phantom 
is  slightly  slower  than  that  predicted  by  the  simulations. 

ers.  Most  manufacturers  are  actively  pursuing  the  development  of 
diagnostic  transducers  capable  of  generating  the  required  intensi¬ 
ties.  While  further  research  is  required  to  investigate  the  thermal 
and  power  issues.  Remote  Palpation  as  proposed  herein  appears 
to  have  clinical  potential. 
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Abstract 

A  method  of  acoustic  Remote  Palpation,  capable  of  imaging  local  variations 
in  the  mechanical  properties  of  tissue,  is  under  investigation.  In  this  method, 
focused  ultrasound  is  used  to  apply  localized  (on  the  order  of  2  mm3)  radi¬ 
ation  force  within  tissue,  and  the  resulting  tissue  displacements  are  mapped 
using  ultrasonic  correlation  based  methods.  The  tissue  displacements  are  in¬ 
versely  proportional  to  the  stiffness  of  the  tissue,  and  thus  a  stiffer  region  of 
tissue  exhibits  smaller  displacements  than  a  more  compliant  region.  A  small 
volume  of  tissue  that  is  significantly  stiffer  than  the  surrounding  medium  (i.e. 
a  lesion)  tends  to  exhibit  rigid  body  motion.  In  this  paper,  the  feasibility  of 
Remote  Palpation  is  demonstrated  experimentally  using  breast  tissue  phan¬ 
toms  with  spherical  lesion  inclusions,  and  in  vitro  liver  samples.  A  single 
diagnostic  transducer  and  modified  ultrasonic  imaging  system  are  used  to 
perform  Remote  Palpation.  The  displacement  images  are  directly  correlated 
to  local  variations  in  tissue  stiffness  with  higher  contrast  than  the  correspond¬ 
ing  B-mode  images  of  the  phantoms.  Relationships  between  acoustic  beam 
parameters,  lesion  characteristics  and  radiation  force  induced  tissue  displace¬ 
ment  patterns  are  investigated  and  discussed.  The  results  show  promise  for 
the  clinical  implementation  of  Remote  Palpation. 

I.  INTRODUCTION 

A.  The  Remote  Palpation  Method 

It  is  hypothesized  that  acoustic  radiation  force  can  be  used  to  generate  localized  tis¬ 
sue  displacements,  and  that  these  displacements  will  be  directly  correlated  with  localized 
variations  in  tissue  stiffness.  It  is  further  hypothesized  that  this  can  be  accomplished  us¬ 
ing  a  single  transducer  on  a  diagnostic  ultrasound  scanner  to  both  generate  the  radiation 
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force  and  detect  the  resulting  displacements.  These  hypotheses  form  the  basis  for  a  new 
imaging  method  called  Remote  Palpation.  In  this  method,  acoustic  radiation  force  is  used 
to  generate  localized  displacements  in  tissue.  These  displacements  are  measured  using  ul¬ 
trasonic  correlation  based  methods.  The  magnitude  of  the  tissue  displacement  is  inversely 
proportional  to  the  local  stiffness  of  the  tissue.  A  volume  of  tissue  that  is  significantly 
stiffer  than  the  surrounding  tissue  (i.e.  a  breast  lesion)  tends  toward  rigid  body  motion, 
and  thus  distributes  the  localized  force  throughout  the  surrounding  tissue.  The  volume  of 
displaced  tissue  is  larger,  and  the  maximum  displacement  smaller  than  that  for  normal  tis¬ 
sue.  Radiation  force  induced  displacement  maps  are  generated  at  multiple  locations,  and 
combined  to  form  a  single  image  of  variations  in  tissue  stiffness  throughout  an  extended 
Field  of  View  (FOV).  A  single  transducer  on  a  diagnostic  scanner  is  used  to  both  generate 
the  high  intensity  ’pushing’  beams  and  track  the  resulting  tissue  displacements. 

B.  Purpose 

Two  potential  clinical  applications  for  Remote  Palpation  are  lesion  detection  and  charac¬ 
terization,  and  the  identification  and  characterization  of  atherosclerosis.  This  paper  focuses 
on  the  former,  as  it  pertains  to  the  early  detection  of  breast  cancer,  which  has  been  shown  to 
significantly  improve  patient  survival.  Existing  methods  of  breast  cancer  detection  include 
screening  mammography  and  palpation,  either  by  patient  self-examination  or  clinical  breast 
exam.  Palpation  relies  on  the  manual  detection  of  differences  in  tissue  stiffness  between 
breast  lesions  and  normal  breast  tissue.  The  success  of  palpation  is  due  to  the  fact  that  the 
elastic  modulus  (or  Young’s  modulus)  of  breast  lesions  is  often  an  order  of  magnitude  greater 
than  that  of  normal  breast  tissue1,2,  i.e.  breast  lesions  feel  ‘hard’  or  ‘lumpy’  as  compared  to 
normal  breast  tissue. 

Differences  in  Young’s  moduli  are  the  basis  for  the  investigation  of  imaging  modalities 
that  provide  information  about  the  stiffness  of  tissue.  Traditionally,  these  have  fallen  into 
two  categories:  1)  Sonoelasticity,  in  which  low  frequency  shear  wave  propagation  is  imaged 
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using  Doppler  or  Magnetic  Resonance  methods.  Estimates  of  the  Elastic  modulus  (Young’s 
modulus)  of  the  tissue  are  based  upon  this  information3-7.  2)  Elastography,  in  which  local 
variations  in  tissue  strain  are  determined  by  measuring  local  displacements  that  occur  during 
global  tissue  compression.  Estimates  of  Elastic  Modulus  are  based  upon  reconstruction 
methods  to  determine  the  elastic  moduli  associated  with  the  measured  strain  fields.  1,8-11 . 
Remote  Palpation  is  similar  to  elastography,  which  has  recently  demonstrated  some  success 
in  the  detection  of  malignant  breast  lesions12,  however  it  has  several  potential  advantages. 
These  include  the  very  localized  application  of  radiation  force  (as  opposed  to  global  external 
compression),  the  decrease  in  maximum  tissue  strain  required  for  lesion  visualization,  and  the 
potential  for  real-time  implementation  without  the  need  for  external  compression  fixtures. 

II.  BACKGROUND 

A.  Acoustic  Radiation  Force 

Acoustic  radiation  force  is  a  unidirectional  force  that  is  applied  to  absorbing  or  reflecting 
targets  in  the  propagation  path  of  an  acoustic  wave.  This  phenomenon  is  caused  by  a  transfer 
of  momentum  from  the  acoustic  wave  to  the  propagation  medium.  The  contribution  of 
absorption  is  in  the  direction  of  wave  propagation,  whereas  the  contribution  of  scattering  is 
dependent  upon  the  angular  scattering  properties  of  the  target.  In  the  event  that  the  target 
is  both  absorbing  and  has  a  scatterer  that  acts  as  a  reflector  with  an  axis  of  symmetry 
perpendicular  to  the  direction  of  wave  propagation  ( e.g .  the  reflecting  target  is  spherical), 
the  radiation  force  is  entirely  in  the  direction  of  wave  propagation.  In  this  situation,  the 
radiation  force  is  given  by13: 

F  =  (no  +  11  s  —  J  7  cos  OrdrdO^  (E) ,  (1) 

where  IIa  is  the  total  power  absorbed  by  the  target,  IIS  is  the  total  power  scattered  by  the 
target,  7  is  the  magnitude  of  the  scattered  intensity,  0  is  the  scattering  angle,  rdrdO  is  an 
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area  element  of  the  projection  of  the  target  onto  the  axial/lateral  plane,  and  (E)  is  the 
temporal  average  energy  density  of  the  propagating  acoustic  wave. 

If  there  is  no  absorption,  and  the  target  is  perfectly  reflecting  ( i.e .  6  —  180°),  the 
integral  term  becomes  — IIS  and  the  total  radiation  force  is  {2IIS  ( E )}.  The  factor  of  2  can 
be  intuitively  explained  by  the  fact  that  the  target  stops  the  wave  from  propagating,  and 
reflects  it  in  the  opposite  direction14.  If  the  target  is  a  Rayleigh  scatterer  (i.e.  its  diameter  is 
much  smaller  than  a  wavelength),  scattering  occurs  in  all  directions  equally,  and  the  integral 
term  is  zero.  Thus  the  total  radiation  force  on  the  scatterer  is  {(IIa  +  IIS)  ( E )}. 

For  a  tissue  model  comprised  of  a  collection  of  Rayleigh  scatterers,  one  can  sum  the 
contribution  of  the  radiation  force  from  each  scatterer  to  determine  the  total  force  due  to 
scattering.  However,  in  tissue,  the  majority  of  the  attenuation  of  an  acoustic  wave  is  due  to 
absorption15,  thus  the  contribution  to  the  radiation  force  by  scattering  in  soft  tissue  can,  in 
general,  be  neglected.  Eq.  1  then  becomes: 


F  =  Ua  (E). 


(2) 


Assuming  plane  wave  propagation,  the  radiation  force  applied  to  tissue  is  thus  14,16  18 


F  = 


2 al_ 
c 


(3) 


where  F  (which  is  in  the  form  of  a  body  force,  or  force  per  unit  volume)  is  the  acoustic 
radiation  force  [kg/(s2cm2)],  or  [dynes/(1000  cm3)],  c  [m/s]  is  the  speed  of  sound  in  the 
medium,  a  [m_1]  is  the  absorption  coefficient  of  the  tissue,  and  I  [Watts/cm2]  is  the  temporal 
average  intensity  of  the  acoustic  beam  at  a  given  point  in  the  tissue.  For  a  focused  acoustic 
beam  propagating  through  soft  tissue,  the  ’target’  is  the  tissue,  and  the  force  is  applied 
throughout  the  focal  region  of  the  acoustic  beam. 

An  important  phenomenon  that  is  not  modeled  by  Equation  3  is  the  enhancement  of 
radiation  force  generated  by  nonlinear  propagation  of  an  acoustic  wave19,18,20,21 .  For  the  same 
temporal  average  intensity,  a  wave  with  higher  pressure  amplitude  and  shorter  pulse  duration 
generates  a  larger  radiation  force  than  does  a  lower  amplitude,  longer  duration  wave.  This 
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is  due  to  the  higher  order  harmonics  generated  by  nonlinear  propagation,  which  result  in  an 
increase  in  absorption19,18,20,21.  Nonlinear  propagation  also  narrows  the  transmit  beam  and 
shifts  the  intensity  field  of  a  focused  acoustic  beam20.  Increases  in  radiation  force  resulting 
from  nonlinear  enhancement  by  a  factor  of  2.6  have  been  observed  in  breast  applications  in 
vivo22.  These  effects  can  be  modeled  using  Equation  4,  where  the  contribution  to  radiation 
force  of  each  harmonic  component  (n)  of  the  acoustic  wave  can  be  computed  at  every  point 
in  the  focal  region  of  the  acoustic  beam  21,18: 

F  =  -J>nIn.  (4) 

c  n 

The  potential  applications  of  acoustic  radiation  force  in  diagnostic  imaging  have  hereto¬ 
fore  been  widely  ignored,  primarily  because  the  available  energy  from  diagnostic  ultrasound 
systems  was  low.  However,  advances  in  transducer  and  system  design  in  recent  years  have 
resulted  in  considerable  increases  in  the  maximum  possible  acoustic  energy  output  from 
these  systems.  This  has  sparked  a  renewed  interest  by  several  laboratories  in  the  potential 
applications  of  acoustic  radiation  force.  Our  laboratory  demonstrated  the  use  of  acoustic 
radiation  force  to  induce  acoustic  streaming  in  cyst  fluid,  and  thus  differentiate  fluid-filled 
from  solid  breast  lesions  in  vivo22,23.  Several  authors  have  proposed  the  use  of  acoustic 
radiation  force  to  remotely  characterize  the  mechanical  properties  of  tissue  19,24~27.  One 
application  involves  the  use  of  a  radiation  force  field  oscillating  at  the  beat  frequency  of 
two  confocal  transducers  to  vibrate  an  object;  the  vibrations  are  detected  by  a  hydrophone, 
and  are  used  to  generate  an  image27.  In  another  application,  acoustic  radiation  force  is 
used  to  displace  tissue,  and  the  speed  of  the  shear  waves  generated  immediately  after  force 
removal  is  monitored  to  characterize  variations  in  tissue  Young’s  modulus19.  In  yet  another 
application,  acoustic  radiation  force  is  used  to  manipulate  the  vitreous  humor  of  the  eye24. 

B.  Tissue  Response  to  Radiation  Force 

Remote  Palpation  applies  localized  forces  to  tissue,  which  can  be  described  for  idealized 
cases  by  elasticity  theory.  For  an  infinitely  small  distributed  force  volume,  the  strain  field 
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can  be  derived  from  the  analytic  solution  for  a  point  load  in  an  infinite  elastic  solid28: 

=  ~(Pz  (1  +  l/)  (4*2  (~l  +  ^)  +  ^2  (-1  +  4*/))) 

8  Etv  (r2  +  z 2)2  (—1  +  v ) 

_  -  (Pz  (— 2r2  +  ,z2)  (!  +  »/))  _  -(Pz  (1  +  t/)) 

8i?7r  (r2  +  z2)2  (—l  +  i/)  8Eiv(r2  +  z2)2  (— 1  +  v) 


where  eZtTye  are  the  3  dimensional  normal  strains.  The  variable  2  is  the  axial  distance  along 
the  line-of-action  of  the  applied  force,  P  is  the  magnitude  of  the  applied  force,  r  is  the  radial 
distance  from  the  point  of  applied  force,  E  is  the  Young’s  modulus  of  the  material,  and  v 
is  Poisson’s  ratio.  The  displacement  of  the  material  in  the  direction  of  the  applied  force  is 
obtained  by  integration  of  e228: 

P  (1  +  v)  (4 *2  (-1  +  u)  +  r2  (-3  +  4  „)) 

w  = - 3 -  (6) 

8  E  tv  (r2  +  z2)2  (—l  +  i/) 

Solving  these  equations  for  a  magnitude  of  force  similar  to  that  used  in  RP  ( i.e .  P  =  0.01 
dynes,  v  —  0.49,  E  =  30,000  dynes/cm2  (3.0  kPa))  gives  an  estimate  of  the  size  of  the 
associated  displacement  and  strain  fields.  The  solution  indicates  that  the  spatial  extent  of 
the  axial  strain  field  is  very  small,  with  strains  decreasing  to  a  negligible  amount  at  a  distance 
of  0.2  mm  from  the  loading  point.  Strains  in  the  radial  and  circumferential  directions  are 
even  smaller.  This  implies  that  lesions  that  are  separated  by  more  than  0.4  mm  can  be 
differentiated  when  excited  with  point  loads  at  a  0.2  mm  spacing. 

Several  insights  can  be  gained  from  this  simplified  analysis.  First,  RP  can  locally  sample 
tissue  stiffness  with  minimal  confounding  artifacts  from  neighboring  tissue  stiffness  variations 
{i.e.  two  lesions  in  close  proximity  will  not  affect  each  other).  Second,  given  that  the  strain 
and  displacement  fields  are  highly  localized,  displacements  measured  at  the  point  of  force 
application  will  be  highly  correlated  with  tissue  stiffness  at  that  location.  Thus,  for  a  local 
increase  in  tissue  stiffness,  smaller  displacements  will  occur  than  in  the  surrounding  tissue. 
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C.  High  Intensity  Acoustic  Beams 

Potential  risks  associated  with  Remote  Palpation  arise  from  the  fact  that  high  intensity 
acoustic  pulses  are  utilized.  The  FDA  provides  two  indices  to  monitor  the  safety  of  diagnostic 
ultrasound.  The  Mechanical  Index  (MI)  is  indicative  of  the  potential  for  cavitation,  and  the 
Thermal  Index  (TI)  is  indicative  of  the  potential  temperature  rise.  The  high  intensity 
pulses  used  in  Remote  Palpation  will  not  exceed  the  MI  limit  of  1.929,  thus  cavitation  is 
not  anticipated  (primarily  because  of  the  relatively  high  frequencies  (7  to  14  MHz)).  Nor  is 
it  anticipated  that  the  TI  will  exceed  6,  which  is  the  limit  above  which  the  FDA  requires 
special  provisions  for  approval  of  use.  Finite  element  models  of  the  heating  associated  with 
the  high  intensity  beams  used  for  Remote  Palpation  indicate  that  the  temperature  rise  will 
be  less  than  1°C30,  which  would  correspond  to  a  TI  of  1. 

Throughout  this  paper,  the  derated  (or  in  situ  estimated)  spatial  peak  temporal  average 
intensity  for  the  different  acoustic  beams  is  provided  in  order  to  calibrate  the  reader.  The 
FDA  currently  limits  the  spatial  peak  temporal  average  intensity  (Ispta.3)  to  0.72  W/cm2  in 
situ29,  which  is  intended  as  an  indirect  indicator  of  potential  tissue  heating.  However,  this 
limit  was  determined  assuming  an  indefinite  application  time.  Short  duration,  high  intensity 
acoustic  pulses  as  are  used  in  Remote  Palpation  (i.e.  10  milliseconds)  were  not  foreseen  in 
the  development  of  this  FDA  limit,  and  thus  the  TI  and/or  temperature  increase  estimates 
provide  a  more  appropriate  indication  of  the  potential  for  thermal  effects. 

D.  Goals 

The  goals  of  the  work  presented  herein  are:  first,  to  experimentally  demonstrate  the 
feasibility  of  Remote  Palpation  imaging;  second,  to  determine  whether  the  displacement 
images  generated  during  Remote  Palpation  imaging  are  directly  correlated  with  variations 
in  tissue  stiffness;  and  third,  to  determine  whether  a  single  transducer  on  a  diagnostic  scanner 
can  be  used  to  perform  Remote  Palpation  imaging. 
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III.  METHODS 

Experiments  were  performed  with  a  Siemens  Elegra  scanner  (Siemens  Medical  Systems, 
Ultrasound  Group,  Issaquah,  WA),  that  has  been  modified  to  allow  user  control  of  the  acous¬ 
tic  beam  sequences  and  intensities,  as  well  as  providing  access  to  the  raw  Radio-Frequency 
(RF)  data.  A  Siemens  75L40  linear  array  was  used  for  these  experiments.  Figure  2  provides 
a  schematic  of  a  typical  linear  array,  and  defines  the  different  dimensions  discussed  below. 
The  75L40  array  consists  of  194  elements,  each  of  which  are  5  mm  tall  and  approximately 
0.2  mm  wide,  with  a  center  frequency  of  7.2  MHz.  The  number  of  active  elements  can  be 
selected  electronically  to  adjust  the  lateral  transmit  aperture  width.  There  is  a  fixed-focus 
acoustic  lens  on  each  element  that  focuses  in  the  elevation  dimension.  Lateral  focusing  is  ac¬ 
complished  electronically  by  applying  the  appropriate  delays  to  each  active  element  (Fig.  2). 
Interrogation  of  different  spatial  locations  is  accomplished  by  using  different  sub-apertures 
(or  groups  of  elements)  within  the  array.  In  all  of  these  experiments,  the  lateral  aperture  was 
considerably  larger  than  the  elevation  aperture.  This  resulted  in  an  asymmetric  focal  region 
that  was  more  tightly  focused  in  the  lateral  dimension  than  in  the  elevation  dimension  (in 
contrast  to  a  piston,  which  has  an  axi-symmetric  focal  region). 

The  tracking  beams  were  standard  diagnostic  B-mode  pulses  ( i.e .  single  cycle  pulses, 
F/l  focal  configuration,  apodized,  Ispta  <  0.1  W/cm2,MI  ~  0.4).  The  pulse  length,  transmit 
voltage,  and  transmit  F-number  (i.e.  number  of  active  transmit  elements)  of  the  beams 
were  varied  in  the  experiments.  The  pulse  lengths  ranged  from  0.8  to  5.8  /rsec,  the  transmit 
voltages  ranged  from  30  to  140  volts  (in  contrast  to  20  volts  for  the  tracking  beams),  the 
F-numbers  ranged  from  F/l  to  F/3,  the  Pulse  Repetition  Frequencies  (PRFs)  were  varied 
from  3  to  10  kHz,  and  they  were  not  apodized. 

Several  beam  sequences  were  designed  to  investigate  the  temporal  and  spatial  response 
of  the  phantoms.  The  experiments  are  divided  into  three  groups:  single  pushing  location 
with  2-D  tracking  region,  single  A-line  interrogation,  and  multiple  pushing  locations. 
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A.  Single  pushing  location  with  2-D  tracking  region 

For  these  experiments,  the  transducer  was  held  in  a  stationary  location.  The  pulse 
sequence  was  designed  to  transmit  a  series  of  13  spatially  distributed  tracking  beams  (i.e. 
B-mode  interrogation  using  13  laterally  spaced  lines,  Fig.  3,  top)  at  a  PRF  of  3584  Hz.  This 
was  followed  by  10  milliseconds  of  pushing  beams  fired  at  a  PRF  of  3584  Hz  along  a  single 
line  of  flight  in  the  center  of  the  tracking  beams,  followed  by  another  series  of  tracking  beams 
interspersed  with  pushing  beams  (every  other  beam,  see  Fig.  3).  The  spacing  of  the  tracking 
beams  was  0.7  mm.  This  sequence  was  repeated  for  50  milliseconds.  The  raw  RF  data  was 
stored  for  off-line  processing. 


B.  Single  A-line  interrogation 

This  sequence  was  designed  to  fire  all  of  the  tracking  beams  in  the  same  spatial  location 
as  the  pushing  beams  (i.e.  A-line  interrogation,  Fig.  3,  bottom).  The  Pulse  Repetition 
Frequency  (PRF)  for  these  experiments  was  10240  Hz.  Every  other  beam  was  a  tracking 
beam,  thus  the  PRF  for  the  tracking  beams  was  5120  Hz.  The  PRF  of  the  pushing  beams 
was  varied  by  firing  either  every  pushing  pulse  (PRF  5120  Hz,  as  shown  in  Fig.  3,  bottom), 
every  fourth  pushing  pulse  (PRF  1280  Hz),  every  fifth  pushing  pulse  (PRF  1024  Hz),  or 
every  sixth  pushing  pulse  (PRF  853  Hz).  The  first  beam  fired  was  always  a  tracking  beam 
which  served  as  a  reference  for  the  initial  position.  The  raw,  Radio-Frequency  (RF)  data 
was  stored  and  processed  off-line. 

C.  Multiple  pushing  locations  with  2-D  tracking  regions 

These  sequences  were  the  same  as  those  described  in  Section  III  A,  however  only  11 
tracking  lines  were  utilized  for  each  pushing  location  (instead  of  13),  and  these  lines  were 
spaced  at  0.09  mm.  This  resulted  in  datasets  from  each  pushing  location  covering  +/-  0.5 
mm  laterally.  Multiple  pushing  locations  were  interrogated  by  connecting  the  transducer 
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to  an  automated  translation  stage,  and  translating  the  transducer  to  interrogate  different 
regions  of  tissue  (using  the  same  acoustic  beam  parameters  and  sequences  at  each  location). 
These  experiments  were  performed  on  two  different  phantoms,  each  with  an  8  mm  diameter 
lesion  located  1.0  cm  from  the  surface  of  the  phantom.  The  phantoms  were  submerged  in 
a  liquid  slurry  phantom  material,  that  was  created  by  blending  a  homogeneous  phantom. 
This  was  done  in  order  to  achieve  the  same  absorption  and  attenuation  of  the  beam  for 
each  axial  transducer  position,  without  applying  compression  to  the  phantom  itself.  The 
pushing  point  locations  spanned  11  mm  laterally,  and  9  mm  axially  with  1  mm  spacing  in 
both  dimensions  (Fig.  8  top).  The  raw  RF  data  from  each  pushing  location  was  stored  for 
off-line  processing. 


D.  Data  Processing 

Off-line  data  processing  was  accomplished  by  performing  cross-correlation  between  se¬ 
quentially  acquired  tracking  lines31,9.  Each  tracking  line  was  divided  into  a  series  of  search 
regions,  and  the  location  of  the  peak  in  the  cross  correlation  function  between  a  kernel  in 
the  first  tracking  line  and  the  corresponding  search  region  in  the  next  tracking  line  was 
used  to  estimate  the  tissue  displacement  in  that  region.  One-wavelength  translations  were 
removed  and  all  displacements  associated  with  a  correlation  coefficient  greater  than  0.99 
were  considered  valid. 


E.  Intensity  measurements 

Pressure  and  intensity  calibration  measurements  were  made  in  accordance  with  the  guide¬ 
lines  provided  by  the  American  Institute  of  Ultrasound  in  Medicine  32 .  These  measurements 
are  complicated  at  higher  transmit  voltages  due  to  acoustic  saturation33.  Saturation  occurs 
at  lower  pressures  in  water  than  it  does  in  phantoms  and  tissue,  thus  accurate  measurement 
of  the  intensities  and  pressures  used  for  the  highest  transmit  voltages  was  not  possible. 
However,  the  relationship  between  acoustic  radiation  force  and  the  resulting  displacement 
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is  linear26.  Therefore,  one  can  obtain  an  estimate  of  the  magnitude  of  the  radiation  force 
generated  by  a  very  high  intensity  acoustic  beam  by  comparing  the  corresponding  peak  dis¬ 
placement  obtained  in  the  same  phantom  using  a  lower  intensity,  quantifiable  beam.  Where 
Eq.  3  applies,  the  ratio  of  the  forces  is  equal  to  the  ratio  of  the  intensities.  Thus  one  can 
obtain  an  estimate  of  the  intensity  in  the  higher  energy  beam  under  the  assumption  of  linear 
propagation  by  evaluating  the  ratio  of  the  peak  displacements. 

The  above  approach  does  not  account  for  the  enhancement  in  radiation  force  caused 
by  nonlinear  propagation20,19,18,22.  This  results  in  an  overestimate  of  the  intensity  of  the 
beam,  due  to  the  assumption  of  a  single  frequency  and  absorption  coefficient  in  Eq.  3.  Even 
so,  this  method  provides  a  reasonable  approximation  of  the  intensities,  and  thus  is  utilized 
to  estimate  the  intensities  used  in  the  experiments  where  actual  measurements  were  not 
possible  ( i.e .  transmit  voltages  of  80  volts  and  higher). 

F.  Phantom  Construction 

Experiments  were  performed  both  in  elastography  phantoms,  and  in  thawed,  de-veined 
calf  liver.  The  elastography  phantoms  were  fabricated  from  gelatin,  graphite,  alcohol,  water, 
and  glutaraldehyde34.  The  lesions  in  the  phantoms  were  all  generated  from  a  single  batch, 
with  a  recipe  corresponding  to  a  higher  Young’s  Modulus  than  the  background  material. 
Table  I  provides  the  recipes  for  the  phantoms.  The  calf  liver  was  purchased  from  the  local 
grocery  store,  frozen,  in  a  vacuum-sealed  package.  It  was  thawed  under  water,  and  the 
experiments  were  performed  on  the  same  day. 

G.  Phantom  Mechanical  Property  Characterization 

Characterization  of  the  Young’s  modulus  of  tissue  and  tissue-like  media  is  very  challeng¬ 
ing.  Only  a  few  reports  exist  in  the  literature  on  the  subject35,36,2.  The  reported  values  vary 
considerably,  possibly  because  Krouskop  et.  al.  applied  pre-compression  while  making  their 
measurements,  whereas  the  other  groups  did  not35,36,2,  and  the  stiffness  of  tissue  increases 
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with  compression2.  The  numbers  provided  here  were  obtained  in  a  method  consistent  with 
those  reported.  This  involved  modeling  of  the  test  apparatus  to  account  for  boundary  con¬ 
dition  deviations  from  the  assumed  theoretical  solution,  and  comparison  of  the  experimental 
measurements  with  both  the  theoretical  model  (which  had  inherent  assumptions)  and  the 
boundary  condition  correction  model34,36.  A  thorough  description  of  these  methods  is  out¬ 
side  of  the  scope  of  this  paper.  It  is  possible  that  our  measurements  suffer  from  a  constant 
offset,  however,  the  relative  stiffnesses  of  the  phantoms  are  consistent  ( i.e .  in  order  from  the 
most  compliant  to  the  most  stiff:  Phantom  C,  Phantom  A,  Phantom  B,  Lesion  Material, 
Phantom  D). 


IV.  RESULTS 

A.  Phantom  Fabrication 

Phantoms  A,  B,  and  C  were  designed  using  the  same  recipe  shown  in  Table  I.  Phantom 
A  had  an  11  mm  diameter  lesion,  Phantoms  B  and  C  each  had  an  8  mm  diameter  lesion. 
Although  the  ingredients  were  the  same  for  each  of  these  phantoms,  the  background  material 
varied  in  its  Young’s  modulus  slightly  between  batches.  This  is  likely  due  to  limited  control 
of  temperature  fluctuations  during  fabrication. 

The  measured  Young’s  modulus  of  the  lesion  material  was  0.5  kPa.  The  Young’s  mod¬ 
ulus  of  the  background  material  in  Phantoms  A,  B  and  C  was  0.07  kPa,  0.13  kPa,  and  .05, 
respectively.  These  values  are  approximately  an  order  of  magnitude  lower  than  those  re¬ 
ported  for  soft  tissue35,36,  however  the  Lesion  to  Tissue  Stiffness  Ratios  (LTSR)  are  between 
3.8  and  10,  which  are  consistent  with  LTSRs  reported  for  breast  tissue35. 

The  Young’s  modulus  of  Phantom  D  was  1.6  kPa.  Phantom  D  was  designed  without  a 
lesion.  The  Young’s  modulus  of  the  liver  sample  used  for  the  experiments  was  not  measured, 
however  reported  values  in  the  literature  range  from  0.4  to  1.7  kPa36.  Comparison  of  the 
steady  state  displacements  generated  with  the  same  Remote  Palpation  imaging  sequence  in 
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both  Phantom  D  and  the  liver  sample  indicate  that  Phantom  D  was  approximately  4  times 
more  stiff  than  the  liver  sample. 

B.  Single  pushing  location  with  2-D  tracking  region 

Fig.  4a  portrays  the  two-dimensional  displacement  profile  generated  in  a  homogeneous 
portion  of  Phantom  A  for  an  F/l  focal  configuration  with  a  lateral  focal  point  of  20  mm, 
and  a  relatively  low  transmit  voltage  (29  volts).  This  image  was  generated  after  10  millisec¬ 
onds  of  force  application,  and  represents  the  tissue  displacement  profile  (in  two  dimensions) 
resulting  from  a  single  force  application.  If  the  acoustic  beam  had  generated  a  point  load, 
one  would  expect  a  small  circular  displacement  profile,  consistent  with  Fig.  1.  However,  the 
focused  acoustic  beam  used  for  these  experiments  is  associated  with  an  intensity  field  that 
is  approximately  the  size  of  its  focal  region.  The  extent  of  the  inner  contour  (-2dB)  shown 
in  Fig.  4a  is  approximately  the  size  of  the  intensity  field  of  the  acoustic  beam  (including 
intensity  values  down  to  -20  dB). 

The  derated  spatial  peak  temporal  average  intensity  of  the  pushing  beams  used  to  gen¬ 
erate  Fig.  4a  was  2.4  W/cm2,  and  the  associated  radiation  force,  as  computed  using  Eq.  3, 
was  100  dynes/cm3.  The  maximum  displacement  is  2.9  microns.  Fig.  4b  shows  the  two- 
dimensional  displacement  profile  generated  when  the  same  beam  sequence  is  focused  in  the 
center  of  an  11  mm  diameter  lesion  in  the  same  phantom.  The  maximum  displacement  is 
smaller  than  in  the  homogeneous  case  (1.4  microns),  and  the  shape  of  the  intensity  field  is 
not  apparent  in  the  contours  of  displacement. 

Fig.  4c  portrays  the  displacement  map  with  the  pushing  beams  focused  in  the  same 
phantom  location  as  in  Fig.  4a,  however  with  a  much  higher  transmit  voltage  (80  volts). 
The  shape  of  the  contours  of  displacement  are  similar  to  Fig.  4a,  however  they  are  shorter 
axially,  and  the  peak  has  shifted  slightly  closer  to  the  transducer.  The  estimated  in  situ  peak 
intensity  (lspta.3)  of  this  beam  is  15  W/cm2  (Section  III  E).  The  maximum  displacement  is 


18  microns. 
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Fig.  4d  was  obtained  using  the  same  beam  sequence  as  in  Fig.  4c,  but  focused  within 
the  lesion  (in  the  same  location  as  in  Fig.  4b).  The  maximum  displacement  is  7  microns. 
Comparison  of  plot  pairs  a-b  and  c-d  indicates  that  their  differences  are  similar.  In  both 
plot  pairs,  the  -10  dB  contours  begin  axially  in  similar  locations.  As  predicted,  the  presence 
of  the  lesion  decreases  the  peak  displacement,  and  alters  the  displacement  profile  from  the 
homogeneous  case.  In  both  Fig.  4b  and  Fig.  4d  the  maximum  displacement  is  approximately 
half  of  that  in  the  corresponding  homogeneous  cases  (Fig.  4a  and  Fig.  4c). 

Fig.  5  portrays  the  displacement  profiles  generated  in  the  calf  liver  sample,  using  different 
focal  positions  (10,  15,  and  20  mm).  In  each  case,  an  F/l  focal  configuration  was  maintained 
(i.e.  the  number  of  active  transmit  elements  was  increased  as  the  focal  position  depth 
increased).  The  estimated  derated  peak  intensity  of  the  beam  used  for  the  20  mm  focus 
is  120  W/cm2  (Section  HIE).  In  each  plot,  the  peak  in  the  displacement  profile  occurs 
slightly  in  front  of  the  focal  position.  Although  the  focal  configuration  was  constant  (F/l), 
the  maximum  displacement  differs  in  each  plot:  12,  11  and  16  microns,  for  plots  a,  b,  and  c 
respectively. 


C.  Single  A-line  interrogation 

Fig.  6  portrays  the  maximum  displacement  through  time  in  Phantom  A.  In  this  sequence, 
all  of  the  tracking  lines  were  fired  along  the  same  line  (spatially)  as  the  pushing  lines,  and 
the  tracking  and  pushing  lines  were  interspersed  (i.e.  track,  push,  track,  push).  The  system 
was  in  an  F/l  focal  configuration,  and  the  transmit  voltage  was  140  volts. 

The  curves  in  the  top  plot  were  generated  using  the  same  PRF  for  both  the  pushing  and 
tracking  pulses  (3584  Hz),  but  different  pulse  lengths.  As  expected,  the  increase  in  steady 
state  displacement  is  linearly  related  to  the  increase  in  pulse  length  (i.e.  the  increase  in 
intensity).  For  example,  the  steady  state  displacement  for  the  3.4  microsecond  pulse  is  double 
that  of  the  1.7  microsecond  pulse.  Also  as  expected,  the  time  constants  for  displacement 
were  similar  for  the  different  pulse  lengths  (~6  milliseconds,  determined  by  fitting  the  curves 
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to  an  exponential). 

The  bottom  plot  was  generated  using  pushing  beams  with  the  same  pulse  length  and 
transmit  voltage,  but  with  different  PRFs  (the  PRF  of  the  tracking  beams  was  held  constant 
at  3584  Hz) .  The  frequency  of  the  ’rectified  sine  wave’  apparent  in  each  of  the  three  curves 
corresponds  to  the  PRF  of  the  high  intensity  beams  for  each  experiment.  This  temporal 
response  is  due  to  the  relaxation  of  the  tissue  between  pushing  beams.  Note  that  the  magni¬ 
tude  of  the  relaxation  increases  with  increasing  mean  tissue  displacement.  This  relaxation  is 
not  apparent  in  the  top  plot  because  the  PRF  of  the  tracking  beams  was  the  same  as  that  for 
the  pushing  beams,  thus  the  relaxation  response  was  undersampled.  Again,  as  expected,  the 
time  constant  for  each  of  these  curves  was  similar  to  those  in  the  top  plot  (~6  milliseconds). 

Fig.  7  portrays  the  temporal  response  of  Phantoms  A,  D,  and  the  liver  sample  during 
force  application  and  after  removal  of  the  force.  The  ramp-down  time  constants  are  4.2,  1.5, 
and  7.9  milliseconds,  respectively.  This  parameter  is  indicative  of  the  damping  present  in  the 
tissue.  The  ramp-up  time  for  Phantom  D  is  about  the  same  (1.6  milliseconds),  however  it  is 
slightly  faster  for  both  Phantom  A  (5.1  milliseconds),  and  the  liver  sample  (6.2  milliseconds). 

D.  Multiple  pushing  locations  with  2-D  tracking  regions 

Fig.  8(top)  provides  a  schematic  showing  the  position  of  the  different  pushing  locations 
(i.e.  focal  points)  used  in  the  multiple  pushing  location  experiments,  and  their  relationship 
to  the  lesion  in  Phantom  C.  The  sub-plots  show  the  displacement  in  the  same  spatial  location 
as  the  pushing  beam  after  10  milliseconds  of  force  application  in  each  grid  location.  Each 
subplot  represents  the  displacements  obtained  at  all  of  the  axial  positions  in  one  lateral 
location  in  the  grid  (as  indicated  by  the  letters  a-f).  The  presence  of  the  lesion  is  clear  in 
these  plots,  where  a  ’dip’  in  the  axial  profile  occurs  in  each  lateral  location  that  intersected 
the  lesion.  The  extent  of  the  ’dip’  corresponds  to  the  length  of  the  interrogated  cross-section 
of  the  lesion. 

Fig.  9a  portrays  the  combination  of  all  of  the  displacement  data  from  the  different  pushing 
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locations  (from  Fig.  8)  into  a  single  displacement  image.  In  this  image,  the  displacement 
at  each  pixel  is  determined  by  computing  the  average  displacement  from  all  of  the  pushing 
locations  whose  focal  regions  overlap  with  that  pixel.  The  lesion  is  clearly  apparent  in  the 
displacement  image,  while  it  is  not  as  evident  in  the  conventional  B-mode  image  (Fig.  9b). 

Fig.  10a  represents  a  combined  displacement  image  of  Phantom  B  that  was  generated  in 
the  same  way  as  Fig.  9a.  Phantom  B  had  a  slight  tear  along  the  interface  between  the  lesion 
and  the  background  on  the  left  side  of  the  image.  This  region  filled  with  the  liquid  slurry 
phantom  material  that  surrounded  the  phantom.  Therefore,  it  was  much  more  compliant 
than  the  phantom  itself.  This  is  clearly  apparent  in  the  displacement  map  shown  in  Fig.  10a, 
as  well  as  in  the  corresponding  B-mode  image  (Fig.  10b). 

V.  DISCUSSION 

Several  insights  can  be  gained  from  the  simplified  point  load  analysis  of  radiation  force 
presented  in  the  background  section  (Fig.  1).  Most  significantly,  the  limited  spatial  extent 
of  the  force  leads  to  two  advantages:  1)  Remote  Palpation  can  locally  sample  tissue  stiffness 
with  minimal  confounding  artifacts  from  neighboring  tissue  stiffness  variations.  This  is  in 
contrast  to  elastography,  where  a  stiff  or  soft  inclusion  impacts  the  strain  field  up  to  2 
diameters  away37.  2)  Given  that  the  strain  and  displacement  fields  are  highly  localized, 
displacements  measured  at  the  point  of  force  application  are  highly  correlated  with  tissue 
stiffness  at  that  location.  For  a  local  increase  in  tissue  stiffness,  smaller  displacements  occur 
than  in  the  surrounding  tissue.  Thus,  the  Remote  Palpation  method  provides  the  ability  to 
generate  high  resolution  maps  of  tissue  displacement  that  are  directly  correlated  with  local 
variations  in  tissue  stiffness  (Figs.  9  and  10). 

Figures  9  and  10  clearly  demonstrate  the  feasibility  of  acoustic  Remote  Palpation  imag¬ 
ing.  The  homogeneous  background  medium  is  a  fairly  uniform  gray  color,  and  the  stiffer 
lesion  is  clearly  apparent  as  a  darker  region.  The  spatial  extent  of  the  lesions  shown  in 
these  images  is  comparable  to  that  in  the  corresponding  B-mode  images.  The  presence  of 
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the  lesion  is  not  apparent  in  the  displacement  maps  beyond  the  lesion  boundaries.  This  is 
clearly  evident  in  Fig.  8,  where  the  displacements  on  either  side  of  the  ’dip’  are  very  similar 
to  those  obtained  in  the  same  axial  location  in  the  absence  of  a  lesion  ( e.g .  subplots  8d 
and  8f).  This  supports  the  hypotheses  that  Remote  Palpation  imaging  can  be  performed 
in  complex  media  without  suffering  from  artifacts  from  neighboring  structures,  and  that 
the  displacements  achieved  in  Remote  Palpation  imaging  are  directly  correlated  to  tissue 
stiffness. 

As  predicted  by  finite  element  simulations26,  the  displacement  profiles  remain  approx¬ 
imately  the  same  shape  for  increases  in  intensity;  and  the  displacement  magnitude  scales 
with  intensity  (Fig.  4a, c).  The  slight  decrease  in  axial  extent  of  the  displacement  profile, 
and  the  slight  shift  of  the  displacement  peak  towards  the  transducer  is  likely  a  result  of  non¬ 
linear  propagation  of  the  higher  pressure  amplitude  acoustic  beam  20,38 .  In  both  cases,  the 
presence  of  the  lesion  results  in  a  considerably  different  displacement  profile,  and  it  appears 
to  move  as  a  rigid  body  in  the  high  intensity  case.  The  ratio  of  the  displacement  at  the  focal 
point  in  the  presence/ absence  of  the  lesion  for  the  low  intensity  case  is  2.3  (Fig.  4a, b);  and 
for  the  high  intensity  case  this  ratio  is  2.6  (Fig.  4c, d).  Previously  presented  finite  element 
simulations  suggest  that  these  numbers  should  be  the  same  26 ,  however,  these  simulations 
did  not  account  for  nonlinear  propagation.  It  is  possible  that  the  lesion  and  the  background 
medium  have  slight  differences  in  nonlinearity,  which  could  account  for  the  slightly  larger 
ratio  in  the  high  intensity  case. 

In  comparing  Figs.  4a  and  c,  and  Figs.  4b  and  d,  it  is  clear  that  the  images  generated  at 
lower  intensities  exhibit  more  noise.  This  is  due  to  the  lower  displacements  in  these  images. 
The  minimum  detectable  displacement  for  a  correlation  based  tracking  algorithm  is  depen¬ 
dent  upon  several  factors:  system  Signal-to-Noise  Ratio  (SNR),  pulse-to-pulse  correlation, 
signal  bandwidth,  as  well  as  the  center  frequency  and  pulse  length  of  the  tracking  beams. 
Utilizing  values  for  the  above  parameters  that  are  applicable  to  the  experiments  performed 
herein  ( i.e .  40  dB,  0.99,  70%,  7.2  MHz,  0.3/xsec,  respectively),  the  lower  limit  of  displacement 
that  can  be  tracked  is  computed  to  be  to  be  0.5pm39.  The  majority  of  the  displacements 
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in  the  low  intensity  cases  (Figs.  4a  and  b)  are  close  to  this  minimum  detectable  value,  thus 
these  displacement  images  have  a  lower  Signal-to-Noise  ratio  (SNR)  than  the  higher  inten¬ 
sity  images  (Figs.  4c  and  d).  Clearly,  larger  displacements  are  preferable  to  achieve  better 
SNRs. 

Figure  5  represents  the  first  in  vitro  Remote  Palpation  displacement  maps  generated  in 
a  tissue  sample.  The  varying  lateral  focal  point  experiments  were  performed  in  order  to 
investigate  the  limitations  presented  by  a  linear  array  with  a  fixed  elevation  focus.  Each 
of  the  three  maps  was  generated  using  an  F/l  focal  configuration  and  the  same  transmit 
voltage  and  pulse  length.  If  the  arrays  were  axi-symmetric  ( i.e .  pistons),  one  would  expect 
identically  shaped  displacement  profiles,  with  decreases  in  peak  displacement  at  deeper  focal 
positions  due  to  the  attenuation  of  the  overlaying  tissue.  However,  because  a  linear  array 
(asymmetric,  5  mm  x  20  mm)  was  used,  the  peak  displacement  actually  occurs  for  the 
deepest  focal  position  (Fig.  5c).  This  focal  position  is  closer  to  the  elevational  focal  point 
than  the  other  two  focal  positions,  which  results  in  increased  intensities  in  this  location. 
Remote  Palpation  will  ideally  be  implemented  by  holding  the  transducer  in  a  stationary 
position,  and  using  electronic  focusing  to  apply  radiation  force  to  multiple  pushing  locations 
within  an  extended  Field  Of  View  (FOV).  The  implications  of  Figure  5  are  that  some  type  of 
depth/focal  position  dependent  scaling  may  be  necessary  to  achieve  a  uniform  displacement 
map  in  a  homogeneous  medium.  This  would  be  similar  to  Time  Gain  Control  (TGC),  which 
is  available  on  diagnostic  scanners  to  account  for  tissue  attenuation  and  focal  gain.  These 
results  also  suggest  that  a  two-dimensional  array  capable  of  electronically  focusing  in  both 
the  lateral  and  elevation  dimensions  is  preferable  to  a  linear  array. 

Evaluation  of  the  steady  state  displacements  shown  in  the  top  of  Fig.  6  indicates  a  linear 
relationship  between  steady  state  displacement  and  acoustic  pulse  length.  Because  the  same 
PRF  and  excitation  voltage  were  used  for  each  pulse,  any  nonlinearities  that  were  present 
were  the  same  for  all  pulse  lengths  (i.e.  each  pulse  had  the  same  peak  pressure  amplitude). 
Thus,  the  ratio  of  the  increase  in  pulse  length,  which  is  approximately  equal  to  the  ratio  of 
the  increase  in  intensity,  is  equal  to  the  increase  in  radiation  force  (Eq.  3).  Therefore,  this 
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plot  portrays  a  linear  relationship  between  radiation  force  and  intensity,  which  was  predicted 
by  finite  element  simulations26. 

The  bottom  plot  in  Fig.  6  was  generated  using  the  same  pushing  pulses,  with  different 
PRFs.  Therefore,  the  intensity  was  lower  for  the  lower  PRFs,  which  resulted  in  smaller 
steady  state  displacements.  Because  the  PRF  of  the  tracking  pulses  was  fixed  at  3584  Hz, 
which  is  above  Nyquist  for  each  of  the  pushing  pulse  PRFs,  the  bottom  plot  in  Fig.  6  portrays 
the  relaxation  of  the  phantom  after  each  individual  pushing  pulse.  This  is  apparent  in  the 
’rectified  sine  wave’  that  appears  on  each  curve.  The  low  points  represent  the  last  track 
beam  fired  prior  to  firing  a  pushing  beam.  The  PRF  of  the  pushing  pulses  is  identical  to  the 
frequency  of  the  ’rectified  sine  wave’  for  each  curve.  It  is  interesting  to  note  that  for  each 
curve,  the  magnitude  of  the  relaxation  increases  with  increasing  mean  displacement  ( i.e .  the 
amplitude  of  the  rectified  sine  wave  increases  with  time) .  In  addition,  the  magnitude  of  the 
’rectified  sine  wave’  is  slightly  larger  for  the  lower  PRFs.  This  is  due  to  the  relaxation  of  the 
tissue  between  pushing  pulses.  The  longer  the  time  between  pushing  pulses,  the  more  the 
tissue  relaxes,  resulting  in  a  net  ’loss’  in  maximum  displacement.  This  suggests  that  a  more 
efficient  method  to  implement  Remote  Palpation  will  be  to  use  a  single,  very  long  pushing 
pulse,  instead  of  a  series  of  shorter  duty  cycle  pushing  pulses  (1%  to  3%  duty  cycles)  as  were 
used  in  these  experiments. 

This  point  is  further  demonstrated  by  evaluating  the  displacement  after  the  same  number 
of  pulses  have  been  fired  at  two  different  PRFs  (Fig.  6,  bottom).  For  example,  at  a  PRF 
of  1280  Hz,  after  6  milliseconds,  8  pushing  pulses  have  been  fired,  and  the  displacement  is 
45  microns.  Correspondingly,  8  pushing  pulses  have  been  fired  at  a  PRF  of  853  Hz  after 
8.5  milliseconds,  and  the  displacement  is  36  microns.  If  one  computes  the  temporal  average 
intensity  used  in  both  scenarios  over  8.5  milliseconds,  it  is  the  same.  However,  by  applying 
the  energy  initially,  instead  of  evenly  spacing  it  over  the  entire  time,  an  increase  of  20%  is 
achieved  in  the  maximum  displacement. 

Figure  7  allows  the  evaluation  of  the  relaxation  time  of  two  different  phantoms  and 
the  liver  sample.  This  quantity  is  indicative  of  the  damping  present  in  the  tissue.  The  liver 
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sample  is  clearly  more  damped  than  either  of  the  phantoms,  as  it  has  a  longer  relaxation  time. 
Phantom  D  exhibits  considerably  less  damping  than  the  other  two,  in  that  its  relaxation 
curve  appears  almost  linear. 

As  stated  in  the  background  section,  Remote  Palpation  is  similar  to  elastography,  how¬ 
ever  it  has  many  potential  advantages.  The  fact  that  the  radiation  force  field  is  localized 
and  can  be  applied  in  selected  remote  locations  is  significant.  Because  the  force  is  applied 
directly  to  the  position  of  interest,  considerably  smaller  forces  are  required  than  in  the  case 
of  global  compression.  Thus  one  can  achieve  lesion  detectability  without  exposing  the  tissue 
to  potentially  damaging  larger  strain  fields.  In  addition,  because  displacement  is  directly 
correlated  with  local  variations  in  stiffness,  the  method  does  not  require  complex  reconstruc¬ 
tion  algorithms.  Finally,  in  contrast  to  Elastography,  the  displacements  generated  during 
Remote  Palpation  are  relatively  small  (microns).  Thus  the  correlation  based  tracking  rou¬ 
tines  are  not  subject  to  problems  associated  with  relative  scatterer  motion  in  the  presence 
of  large  strain37. 

Clinical  implementation  of  Remote  Palpation  will  be  similar  to  conventional  mixed 
modes,  such  as  B-Mode/Color-mode.  The  physician  will  identify  a  smaller  FOV  within 
a  B-mode  image  where  Remote  Palpation  will  be  implemented.  Then,  upon  entering  the 
Remote  Palpation  mode,  the  physician  will  hold  the  transducer  in  a  stationary  position  while 
the  FOV  is  interrogated  via  electronic  focusing  and  steering.  Depending  upon  the  size  of  the 
FOV,  the  spacing  of  the  pushing  locations,  and  the  time  spent  at  each  location,  frame  rates 
as  high  as  2  to  5  frames  per  second  may  be  achieved  while  interrogating  a  1.5  cm  square 
FOV.  This  frame  rate  will  allow  the  superposition  of  Remote  Palpation  displacement  images 
with  the  conventional  B-mode  image  in  semi-real-time,  which  will  aid  clinical  evaluation. 

The  resolution  achievable  by  Remote  Palpation  imaging  systems  will  depend  upon  a 
variety  of  factors,  including:  the  size  and  spacing  of  the  pushing  beams,  the  resolution  and 
spacing  of  the  tracking  beams,  the  target  tissue  mechanical  properties  (Equations  5  and  6), 
and  the  methods  of  image  reconstruction  utilized.  The  impact  of  these  factors  on  resolution 
is  under  investigation,  however  it  is  outside  the  scope  of  this  paper. 
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Given  the  comparably  high  intensities  and  unique  pulsing  sequences  of  Remote  Palpation 
imaging,  the  safety  of  the  method  is  of  concern.  The  experiments  presented  herein  indicate 
that  Remote  Palpation  will  require  in  situ  intensities  of  up  to  300  W/cm2  at  each  pushing 
location.  This  is  well  in  excess  of  the  FDA  limit  (0.72  W/cm2)29.  However,  these  high 
intensity  beams  need  only  be  generated  for  time  scales  on  the  order  of  milliseconds.  In  this 
short  time  period,  the  heating  that  would  be  generated  is  less  than  1°C30,  and  is  limited 
spatially  to  a  region  that  is  slightly  smaller  than  the  focal  region  of  the  acoustic  beam. 
Therefore,  assuming  that  adjacent  pushing  beams  do  not  significantly  overlap,  the  maximum 
expected  temperature  rise  will  be  less  than  1°C.  It  is  widely  accepted  that  temperature 
increases  of  less  than  1°C  do  not  pose  a  danger  to  the  patient40.  Thus,  Remote  Palpation 
imaging  as  discussed  herein  should  not  pose  a  danger  to  the  patient. 

VI.  CONCLUSION 

The  work  presented  herein  clearly  demonstrates  the  feasibility  of  acoustic  Remote  Pal¬ 
pation.  Displacement  maps  that  are  directly  correlated  with  variations  in  tissue  stiffness 
were  generated  and  presented.  It  is  possible  to  perform  Remote  Palpation  using  a  single 
transducer  on  a  diagnostic  scanner  to  both  generate  the  required  high  intensity  beams,  and 
track  the  resulting  displacements.  This  will  facilitate  ease  of  implementation  on  transition 
to  the  clinical  setting.  The  intensities  that  are  required  to  induce  detectable  displacements 
in  vivo  (i.e.  tens  of  microns)  will  be  on  the  order  of  200-300  W/cm2,  and  the  required  dwell 
time  less  than  5  milliseconds.  Due  to  the  short  duration  of  exposure,  these  beams  should  not 
cause  appreciable  tissue  heating.  Given  these  findings,  acoustic  Remote  Palpation  clearly 
exhibits  clinical  potential. 
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TABLES 


Gelatin 

Water 

N-propanol 

Graphite 

Glutaraldehyde 

E 

(g) 

(ml) 

(ml) 

(g) 

(25%  soln,  ml) 

(kPa) 

Tissue  (A,B,C) 

13.0 

230.0 

18.0 

18.75 

52.0 

~0.1 

Lesion 

17.0 

222.0 

18.0 

18.75 

60.0 

0.5 

Tissue  (D) 

20.0 

222.0 

18.0 

18.75 

60.0 

1.6 

TABLE  I. 


Table  I:  Phantom  fabrication  recipes  for  300  milliliters  of  solution  34 
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Figure  1:  Percent  axial  strain  (top)  and  associated  axial  displacement  (bottom)  resulting  from  the 
application  of  a  point  load  with  a  force  magnitude  similar  to  that  used  for  Remote  Pal¬ 
pation  (i.e.P=0.01  dynes),  in  an  elastic  medium  with  material  properties  consistent  with 
those  of  tissue  ( i.e.u=0A9 ,  E=30,000  dynes/cm2  (3.0  kPa))28.  Note  that  the  axial  extent 
of  both  fields  has  decreased  to  an  insignificant  amount  within  0.2  mm  of  the  location  of 
force  application  (0mm). 

Figure  2:  Linear  transducer  configuration.  Note  that  there  is  only  1  row  of  elements,  and  each 
element  is  much  taller  (elevation  dimension)  than  it  is  wide  (lateral  dimension).  There  is 
an  acoustic  lens  that  focuses  in  the  elevation  dimension,  and  electronic  delays  are  used  to 
focus  in  the  lateral  dimension.  Thus,  the  focal  point  can  be  varied  in  the  lateral  dimension 
(as  was  done  in  Fig.  5),  whereas  the  elevation  focal  point  is  fixed. 

Figure  3:  Schematic  of  the  beam  sequences  for  the  different  experiments:  a)  single  pushing  location 
with  2-D  (B-mode)  tracking  region,  b)  single  A-line  interrogation.  The  numbers  show  the 
spatial  relationship  between  the  different  lines,  and  the  sequences  indicate  the  temporal 
order  in  which  the  lines  were  fired.  ’7p’  refers  to  a  pushing  pulse  that  is  fired  in  spatial 
location  7.  The  different  spatial  locations  are  interrogated  by  using  different  sub-apertures 
(or  groups  of  elements)  in  the  array. 
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Figure  4:  Displacement  maps  after  10  milliseconds  of  force  application  in  phantom  A.  For  each 
of  the  images,  the  black  contours  represent  levels  of  -2  dB  and  -10  dB  of  the  peak 
displacement,  the  lateral  focal  point  was  20  mm,  and  the  transducer  was  located  at 
the  top  of  the  images,  a)  Displacement  profile  generated  in  a  homogeneous  region 
of  the  phantom  for  an  F/l  focal  configuration  using  a  single  pushing  location  and 
a  relatively  low  transmit  voltage  (29  volts),  and  a  pulse  length  of  2.2  microseconds, 
b)  Displacement  profile  generated  using  the  same  beam  sequence  as  in  a),  with  the 
pushing  location  centered  in  an  11  mm  diameter  lesion.  The  lesion  boundaries  are 
highlighted  in  white.  Both  a)  and  b)  have  the  same  colorbar  scales,  ranging  from 
0  to  3  microns,  c)  Displacement  profile  generated  in  the  same  location  as  a),  also 
using  an  F/l  focal  configuration  but  with  a  higher  transmit  voltage  (80  volts),  and  a 
pulse  length  of  0.8  microseconds.  Note  that  the  displacement  profile  is  similar  to  a), 
however  it  is  shorter  axially,  with  the  peak  shifted  slightly  closer  to  the  transducer. 
This  is  to  be  expected  in  the  presence  of  nonlinear  propagation,  d)  Displacement 
profile  generated  in  the  same  location  as  b),  using  the  same  beam  sequence  as  in  c). 
Both  c)  and  d)  have  the  same  colorbar  scales,  with  displacements  ranging  from  0  to 
18  microns. 

Figure  5:  Displacement  maps  generated  in  calf  liver  after  10  milliseconds  of  force  application. 

For  each  of  the  images,  the  black  contours  represent  levels  of  -1  dB,  -3  dB  and  -6  dB 
of  the  peak  displacement.  The  units  of  displacement  on  the  colorbar  are  microns. 
The  transducer  was  located  at  the  top  of  the  images.  A  constant  focal  configuration 
of  F/l  was  used  to  generate  these  images,  with  a  transmit  voltage  of  140  volts, 
and  a  pulse  length  of  3.3  microseconds,  a)  lateral  focal  point  at  10  mm.  b)  lateral 
focal  point  at  15  mm.  c)  lateral  focal  point  at  20  mm.  Note  that  the  peak  in  the 
displacement  map  moves  with  the  lateral  focal  point,  and  is  slightly  closer  to  the 
transducer  than  the  focal  point.  This  is  expected  with  nonlinear  propagation. 
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Figure  6:  Displacement  through  time  near  the  position  of  peak  displacement  (17.5  mm  from 
the  transducer)  in  phantom  A.  Top:  The  pulse  length  of  the  pushing  pulses  was 
varied  for  each  curve  as  indicated  in  the  legend  in  units  of  microseconds,  but  the 
PRF  was  held  constant  (3584  Hz),  as  was  the  focal  configuration  (F/l)  and  transmit 
voltage  (140  volts).  Note  that,  as  expected,  the  steady  state  displacement  for  the 
3.4  microsecond  pulse  is  double  that  of  the  1.7  microsecond  pulse.  Bottom:  For 
each  of  these  curves,  the  same  pushing  pulses  were  used  (F/l  focal  configuration, 
a  pulse  length  of  5.8  microseconds,  and  a  transmit  voltage  of  140  Volts)  but  with 
different  pushing  pulse  PRFs.  The  frequency  of  the  ’rectified  sine  wave’  apparent  in 
each  of  the  three  plots  corresponds  to  the  PRF  of  the  high  intensity  beams  for  each 
experiment  (as  indicated  in  the  legend,  in  Hz).  This  is  due  to  the  relaxation  of  the 
tissue  between  pushing  beams.  Note  that  the  magnitude  of  the  relaxation  increases 
over  time,  with  increasing  mean  tissue  displacement. 

Figure  7:  Displacement  over  time  in  the  focal  region  of  the  transducer  for  Phantoms  A  and  D 
in  a  homogeneous  portion  of  each  phantom,  as  well  as  in  the  liver  sample.  In  both 
Phantom  A  and  the  liver  sample,  the  pushing  pulses  were  3.3  microseconds  long, 
with  an  F/l  focal  configuration  and  a  transmit  voltage  of  140  volts.  In  Phantom 
D,  the  pushing  pulses  were  5.8  microseconds  long,  with  an  F/l  focal  configuration 
and  a  transmit  voltage  of  140  volts.  In  all  experiments  the  pushing  pulses  and 
tracking  pulses  were  interspersed,  and  the  pushing  pulses  were  fired  for  a  total  of 
52  lines,  and  then  turned  off.  The  PRF  in  the  liver  sample  was  5120  Hz,  whereas 
in  phantoms  A  and  D  it  was  3584  Hz.  The  time  at  which  the  pushing  pulses  were 
turned  off  has  been  aligned  in  these  plots  for  ease  of  analysis.  The  dashed  lines 
represent  exponential  fits  to  the  data.  The  time  constants  for  the  ramp  up  time  are 
5.1,  1.6,  and  6.2  milliseconds  for  phantoms  A,  D,  and  the  liver  sample,  respectively. 
The  time  constants  for  the  relaxation  are  4.2,  1.5,  and  7.9  milliseconds  for  phantoms 
A,  D,  and  the  liver  sample  respectively. 
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Figure  8:  Top  left:  Grid  of  pushing  locations  (9x11  mm,  1mm  spacing  in  both  the  axial 
and  lateral  directions)  and  their  relationship  to  the  lesion.  Each  dot  at  the  grid 
intersections  represents  a  pushing  location  focal  point.  Although  tracking  data  was 
obtained  throughout  the  box  to  either  side  of  each  pushing  location,  only  the  central 
tracking  line  (which  is  co-located  with  the  pushing  line)  is  shown.  Top  right:  Raw 
data  from  fourth  row  of  pushing  locations  (including  all  of  the  lateral  tracking  lines) 
used  to  generate  the  image  shown  in  Figure  9.  Note  that  although  the  focal  position 
of  the  pushing  beams  is  at  a  fixed  range  (17mm),  the  lesion  is  well  defined  from  12 
to  21  mm  axially.  Bottom:  Each  subplot  represents  the  displacement  along  the  axis 
of  the  high  intensity  pushing  beam  for  the  9  axial  pushing  locations  at  each  lateral 
position  shown  on  the  grid  above  (a-f).  The  lesion  is  clearly  apparent  as  a  ’dip’  in 
the  displacement  profiles  in  lines  a-d,  and  the  width  of  the  ’dip’  is  consistent  with 
the  cross-sectional  dimension  of  the  lesion  at  each  lateral  location. 

Figure  9:  a)  Displacement  image  created  from  multiple  pushing  locations  in  phantom  C,  corre¬ 
sponding  to  the  raw  data  shown  in  Fig.  8.  The  pushing  pulses  were  1.7  microseconds 
long,  with  an  F/l  focal  configuration  and  a  transmit  voltage  of  140  Volts.  Their 
PRF  was  3584  Hz.  The  image  was  created  by  combining  the  displacements  in  each 
spatial  location  after  10  milliseconds  of  force  application.  The  units  of  displacement 
on  the  colorbar  are  microns,  b)  Corresponding  B-mode  image. 

Figure  10:  a)  Displacement  image  created  from  multiple  pushing  locations  in  phantom  B.  This 
phantom  had  a  tear  between  the  lesion  and  the  background  material  on  the  left 
side  of  the  lesion.  This  resulted  in  some  of  the  liquid  phantom  leaking  into  the 
tear,  and  hence  a  much  more  compliant  region  on  in  the  vicinity  of  the  tear.  The 
pushing  pulses  were  0.83  microseconds  long,  with  an  F/l  focal  configuration  and 
a  transmit  voltage  of  140  Volts,  and  a  PRF  of  3584  Hz.  The  force  applied  in  this 
experiment  was  half  of  that  used  to  generate  the  image  in  Fig.  9.  The  image  was 
created  by  combining  the  displacements  in  each  spatial  location  after  10  milliseconds 
of  force  application.  The  units  of  displacement  on  the  colorbar  are  microns,  b) 
Corresponding  B-mode  image. 
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Linear  Transducer  Configuration:  1  line 


FIG.  2. 
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B-mode  Image  Lines: 

iiiiilJ]]  10111213 

2-D  Tracking  Region  Sequence  (13  lines): 

1,2,  3,  4,  5,  6, 7,  8,  9,  10, 11,  12, 13,  (track) 

7p,  7p,  7p,  7p, ...  ,7p,  7p,  (7p  =  push  at  line  7) 

1,  7p,  2, 7p,  3,  7p,  4, 7p,  5,  7p,  6, 7p,  (interspersed  track/push) 
7, 7p,  8, 7p,  9,  7p,  10,  7p,  1 1, 7p,  12,  7p,  13, 

7p,  7p,  7p,  7p, ...  (repeat  push  and  interspersed  track/push) 

Single  A-line  interrogation  (only  line  7): 

7, 7p,  7, 7p,  7,  7p, ... 
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FIG.  4. 
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lateral  position,  mm 

a  b  c 


FIG.  5. 
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FIG.  9. 


Nightingale  et.  a L,  JASA 


